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Abstract 

Using a rigorous quantum model a comprehensive study of physical properties of entangled photon pairs generated in 
spontaneous parametric down-conversion in chirped periodically-poled crystals is presented. Spectral, temporal, as well as 
spatial characteristics of photon pairs are analyzed. Spectral bandwidths, photon-pair flux, and entanglement area can be 
effectively controlled by chirping. Quantification of entanglement between photons in a pair is given. Splitting of entanglement 
area in the transverse plane accompanied by spectral splitting has been revealed. Using the model temperature dependencies 
of the experimental intensity profiles reported in literature have been explained. 



I. INTRODUCTION 

The beginning of nonlinear quantum optics can be 
dated back to the sixties of the 20th century, when the 
first nonUnear optical phenomena were observed [H ow- 
ing to the discovery of the laser. The next milestone was 
reached twenty years later by Hong, Ou and Mandel who 
observed strong quantum correlations between photons 
in a photon pair generated in the process of spontaneous 
parametric frequency down-conversion Q . These photon 
pairs represent the elementary quantum objects that can 
be correlated in time, space, momentum, energy, polar- 
ization or even orbital angular momentum. They have 
been used in many applications including tests of Bell's 
inequalitiesJS, quantum cryptography 01 , quantum tele- 
portation jaS], and metrology Q- 

Many different structures emitting entangled photon 
pairs have been developed. NonHnear bulk crystals Q 
and nonlinear wave-guides belong to standard sources of 
these photon pairs nowadays. Also new sources based 
on photonic-crystal fibers j9|, HOl] , layered photonic struc- 
tures nonhnear cavities and photonic 
nanostructures like quantum dots have been investigated. 

However, periodically-poled crystals are distinguished 
from other photon-pair sources for two reasons [3, 
[isl l . High photon-pair fluxes are reached because of the 
use of materials with high nonlinearities due to quasi- 
phase-matching. On the other hand, chirping allows to 
generate photon pairs with extremely wide spectral band- 
widths and subsequently extremely sharp temporal fea- 
tures. Such photon pairs are useful in many experiments, 
e.g., in quantum optical coherence tomography [l^ or de- 
termination of tunnelling times [131 ■ Wide spectral band- 
widths of entangled photon pairs make them prospective 
for parallel processing of quantum information. 

In theory, description of chirped periodically-poled 
crystals has been restricted to simple more-less analyti- 
cally soluble models up to now [3, [3| ■ Properties 
of photon pairs in the transverse plane have not been 
addressed in these models. In the paper, we investigate 
spectral, temporal as well as spatial properties of the 
emitted photon pairs in detail using a complete theory. 
Special attention has been devoted to the role of chirping 



parameter that is found extraordinarily useful in tailor- 
ing properties of photon pairs. Splitting of entanglement 
area that is accompanied by spectral splitting is predicted 
in this model. Also temperature dependencies of inten- 
sity profiles experimentally observed in have been 
addressed using the model. 

The model of spontaneous parametric down-conversion 
is introduced in Sec. II. Quantities characterizing photon 
pairs are defined in Sec. III. Spectral and temporal prop- 
erties as well as entanglement of photons in a pair are 
discussed in Sec. IV. Properties of photon pairs in the 
transverse plane including entanglement area and tem- 
perature dependencies are analyzed in Sec. V. Conclu- 
sions are drawn in Sec. VI. 



II. MODEL OF SPONTANEOUS PARAMETRIC 
DOWN-CONVERSION 

We use multi-mode description of spontaneous para- 
metric down-conversion [2l|. Each field, i.e. pump, signal 
and idler, is represented by plane-waves with given po- 
larizations. Nonlinear process is described by first-order 
perturbation solution of Schodinger's equation in the in- 
teraction picture. Interaction Hamiltonian Hint is in the 
form: 

iJint (t) = Sin^eo J ^ dzx^^'^ (z) : E(+) {z,t) 

X e(-) (z,t)E|"' (z,t)-f H.c. (1) 

Symbol S denotes transverse area of the pump beam, 
eo is permittivity of vacuum, and x^^'' (^) represents a 
third-order tensor of second-order nonlinear susceptibil- 
ity that is assumed wavelength-independent. Symbol : 
is a shorthand for tensor x'^^ with respect to its three 
indices and H.c. replaces a Hermitian conjugated term. 
Symbol L denotes length of the nonlinear crystal. The 
strong pump field is described classically whereas the 
weak signal and idler fields are treated quantally using 
operator amplitudes. An efficient nonlinear interaction 
is observed provided that x, y, and z components of the 
wave vectors fulfil at least approximately phase match- 
ing conditions in the corresponding directions. Assuming 



plane-wave pumping and sufficiently wide crystal perfect 
phase matching in the transverse plane (spanned by x 
and y variables) occurs. 

Positive frequency part of pump-field amplitude E^-* 
has the form: 
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The abbreviation TM (TE) stands for transversal mag- 
netic (electric) waves, ujp is pump-field angular frequency, 

kp^^ the corresponding wave- vector, and Ept))(z,Wp) 
spectrum of the positive-frequency part of the pump-field 
amplitude. Strong coherent pumping is assumed and so 
the depletion of the pump field can be neglected. 

The signal and idler fields are described by negative- 
frequency parts of electric-field amplitude operators eI ■* 
and E , respectively, that can be written as follows [2^ : 
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X a]^^ (ujj) exp [-i ([kj- 5]^ z - ujjt)] , 



j = s,i. Symbol ej^^ (ujj) stands for polarization vector 
of a mode with frequency Uj and polarization ^ whereas 



(ujj) represents the creation operator of this mode. 
Symbol h denotes reduced Planck's constant, n^{iLjj) is 
the index of refraction and c is speed of fight. 

First-order solution of Schodinger's equation gives the 
wave-function of an entangled two-photon state in the 
output plane of the crystal: 
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where x*^^^ is Fourier transform of function x*-^-* (z) and kn 
denotes a grid vector of poHng. Alternatively, the crystal 
length can be given indirectly by number Nl of layers. 
Normalization function g {ujs,i^i) has the form: 



(7) 



Variations of second-order nonlinear susceptibility x'^^ 
along the z-axis can be described using the formula: 



x'^' (z) = Xo^^sign[cos (k„oz - Cz^)]rect 



(8) 



Symbol fc„o denotes a basic grid vector of poling and 
kno = tt/^o where Zq represents the length of the basic 
layer. Parameter C, characterizes chirping of the poled 
structure. Function rect(a;) equals 1 for > x > 1 and 
is zero otherwise. We note that periodical poling can 
be reached by applying strong static electric fields to do- 
mains (defined by contacts prepared by lithography) that 
invert the sign of nonfinear susceptibility x'^^ , I 7|. An 
other method is based on mapping an optical-field in- 
terference pattern to x'^^ nonlinearity that can be done 
provided that a nonfinear material is placed into a strong 
uv radiation l2al. 



where \^oq) denotes the initial signal- and idler-field vac- 
uum state. Contributions described by higher orders of 
the perturbation solution are assumed to be small due 
to relatively low pumping intensity. Taking into account 
Eqs. (HI), lIHl), and (0]) the wave function jV'*^^-') can 
be obtained in the form: 




X Aap {i^s, i^i) [ujs) a\ p [lji) I-^oo) ■ (5) 

Here, Aaf3 (ujs,uji) denotes a two-photon spectral ampli- 
tude giving probability amplitude of emitting a signal 
photon at frequency and polarization a together with 
an idler photon at frequency tOi and polarization [3. Two- 
photon amplitude Aap can be expressed as [23 . [23 . [25| 



A. Simplified description using uniform poling 

Formation of the spectral shape by interference from 
all layers can be studied analytically for uniform poling 
and cw pumping provided that we restrict ourselves to a 
scalar modes and consider propagation along the z axis 
only. The two-photon amplitude A can then be expressed 
as: 

A{ujs,u}i) = (ujs,^i)E'^\i-^po)5{ujpo -uJs- UJi) 

X / dzx^^^z) exp[iAkz]; (9) 

Ak = kp — kg — ki. Integration along the crystal in Eq. ^ 
can be replaced by the sum of contributions coming from 
each layer and then the two-photon amplitude A takes 
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the form: 
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Symbol lo introduced in Eq. (fTOl denotes the length of 
one layer {lo = L/Nl). Function sine in Eq. (fTO| origi- 
nates in phase mismatch of one typical layer whereas su- 
perposition of down-converted fields from all layers gives 
the remaining part of the formula in Eq. I|10p that resem- 
bles the well-known expression describing interference of 
a finite number of waves. Maximum value of the two- 
photon amplitude A in Eq. (fTOl) is reached if the sin func- 
tion in the denominator is zero. This happens under the 
usual quasi-phase-matching condition Afc = 7r(n + l)/lo, 
n being an integer. The number of emitted photon-pairs 
at the corresponding frequencies is then proportional to 
iV|. The expression in Eq. (fTO|) also indicates that the 
spectral width of the signal or idler field is proportional 
to 1/Nl for cw pumping. 



III. CHARACTERISTICS OF PHOTON PAIRS 
EMITTED IN PERIODICALLY POLED CRYS- 
TALS 

Number of generated photon pairs N (per 1 s in cw 
regime) can be easily calculated from two-photon spectral 
amplitude A given in Eq. ([6]): 
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Energy spectrum Ss,i (t^s,i) for the signal or idler field is 
determined similarly with the formula: 



Ss,i{0Js,i) ^ lVjJs,i I duJt^s\A{uJsTl^i) \ ■ (12) 
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Schmidt's decomposition of the two-photon spectral 
amplitude A into the dual basis can be used to quantify 
entanglement between signal and idler photons. Two- 
photon spectral amplitude A{ujs,'^i) can then be ex- 
pressed as 



A{LOs,i^i) = ^An^/^n {uj s) {uJi) , 



(13) 



where A„ are eigenvalues corresponding to base functions 
'0n and 0„ . The base functions V'n and 0„ are determined 
as solutions of the integral equations: 



Ks (W, Uj) tpn (lu) diU = X^tpn (uj) , (14) 



Ki {uj, lo) (f>n {lo) duj = A.^,0„ {uj) . (15) 



Entropy of entanglement E [23| defined as: 
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is used for quantification of entanglement. The higher 
the entropy E, the stronger the entanglement. The num- 
ber of effective independent modes useful in Schmidt's 
decomposition is given by cooperation parameter K [23 |: 
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Temporal properties of photons in a pair can only 
be measured indirectly using interferometers. Relative 
time delay between photons in a pair can be deter- 
mined in Hong-Ou-Mandel interferometer [2§|. Normal- 
ized coincidence-count rate i?„ (t) in this interferometer 
is expressed in the form: 



Rn (r) = l-p(r) 



(18) 



and 
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(19) 



Re denotes the real part of an argument. Symbol Ro is 
a normalization constant; Rq — limT-^oo PiT)- 



IV. TEMPORAL AND SPECTRAL PROPERTIES 
OF PHOTON PAIRS, ENTANGLEMENT 

The analysis is done for a periodically-poled nonlin- 
ear uniaxial crystal LiNbOa. Having in mind optical- 
fiber applications we consider cw plain-wave pumping 
at central wavelength Apo = 752.5 nm. The sig- 
nal and idler fields have wavelengths Aso = Aio = 
1505 nm (Aso — 1392.1 nm and A^o = 1637.8 nm) 
in the spectrally-degenerate (spectrally-non-degenerate) 
case. Assuming uniform poling the length identical for 
all layers is determined such that quasi-phase matching 
is reached for TM polarizations of all interacting fields 
at given central frequencies assuming collinear propaga- 
tion. Since the studied material is anisotropic and we 
consider the optical axis parallel to the y axis, polariza- 
tions of all interacting fields correspond to extra-ordinary 
waves in collinear geometry. The length of layers 
for the spectrally-degenerate (spectrally-non-degenerate) 
case equals 8.9000 ^m (8.9286 ^m). 

If chirping in poling is assumed, the length of layer 
in the middle of the structure equals that for uniform 
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poling. The value of used chirping parameter C (C = 
1 X 10"^ Min~^) corresponds to the case in which shorter 
layers occur at the beginning of the crystal whereas 
longer layers are at the end of the crystal . We note 
that if the value of chirping parameter ^ is sufficiently 
small, layer lengths depend more-less linearly on their 
positions. We also note that the effect of multiple re- 
flections at crystal boundaries that leads to fast spec- 
tral oscillations (typical for a Fabry-Perot etalon [30|) 
occurs in periodically-poled crystals. However, we sup- 
press them in our analysis by making local averaging of 
the signal- and idler-fleld spectra. This is in accord with 
usual experimental conditions in which these oscillations 
are smoothed out by a final spectral resolution of detec- 
tors. 

Note that, the Fourier transform of Eq. ^ contains 
several spectral components. The component +kn as- 
sures phase-matching at the central frequencies men- 
tioned above. However, the component with — A:„ also 
exists but cannot provide phase matching in the consid- 
ered frequency range for the studied nonlinear crystal and 
so we omit it in further investigations. 



A. Spectrum and number of generated photon 
pairs 
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Typical spectra for both frequency-degenerate and 
frequency-non-degenerate down-conversion for uniform 
as well as chirped poling are shown in Fig. [TJ 

Fast modulation in these spectra has its origin in inter- 
ference of two-photon amplitudes coming from different 
nonlinear layers of the crystal as expressions in Eqs. lfTO|) 
and (fT2|) indicate. In frequency-non-degenerate case and 
for uniform poling, there occur two peaks centered at 
different frequencies (see Fig. [lb) . These contributions 
mutually interfere in the frequency region where they 
overlap. Chirped poling leads to wide broadening of 
the signal- and idler-field spectra, as demonstrated in 
Fig.[T^,b. The reason is that two-photon amplitudes from 
different layers have shifted ranges of the emitted signal- 
and idler-field frequencies and so when two-photon am- 
plitudes from all layers are superposed, a much broader 
range of frequencies is covered. However this broaden- 
ing is at the expense of lowering photon-pair emission 
rates. This occurs because nonlinear layers of different 
lengths and generating different two-photon amplitudes 
do not admit perfect constructive interference at opti- 
mum frequencies (as it occurs for uniform chirping). We 
also note that wider spatial emission angles occur for the 
higher values of chirping and so a considerable portion 
of the emitted signal does not contribute to the studied 
collinear case. 

Successive forming of the signal-field spectrum as the 
number Nl of layers increases is depicted in Fig. for 
uniform poHng. Already two layers (that form the basic 
element of the structure) provide a spectrum that prefers 
the area around the central frequency. As the curves in 




FIG. 1: Energy spectrum Sa for a) spectrally degenerate and 
b) spectrally non-degenerate down-conversion. Curves de- 
noted as U (C) are appropriate for uniform (chirped) poling 
of the crystal; Nl = 1000 in a) and b). Forming of nor- 
malized energy spectrum 5"°"^™ with the increasing number 
Nl of layers is documented in c) for uniform poling, spectral 
degeneration, and Nl = 1,2, 50, and 1000. 



Fig. Hfc document the larger the number Nl of layers 
the narrower the spectrum [compare the expression in 
Eq. idOD]. 

Number N of generated photon pairs and signal- and 
idler-field spectral widths are the most important param- 
eters of a photon-pair source in many appHcations. We 
show their dependence on the number N^ of layers (giv- 
ing length L of the crystal) and chirping parameter C in 
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FIG. 2: (Color online) Contour plots of a) number of gen- 
erated photon pairs and b) spectral width ASs (FWHM) of 
the signal field as functions of the number of layers A''^ and 
chirping parameter (^. Frequency degenerate case is assumed. 



Fig. [21 We can observe the following general tendencies. 
Provided that ^ is fixed increasing number Nl of layers 
gives higher numbers N of generated photon pairs but 
lowers spectral widths of the signal and idler fields. On 
the other hand and keeping Nl fixed increasing values 
of chirping parameter ( lead to decreasing of the num- 
ber N of generated photon pairs and broader signal- and 
idler-field spectra. 

Defined spectral widths of the signal and idler fields 
are required in many applications. The spectral widths 
are determined by the lengths of the first and last nonlin- 
ear layers in a crystal with chirped poling. Number TV of 
generated photon pairs increases naturally with increas- 
ing number Nl of layers. Since lengths of the first and 
last layers are fixed, a larger number Nl of layers means 
a lower value of chirping parameter (. Our numerical 
investigations have revealed that the product of number 
N of generated photon pairs and chirping parameter C is 
roughly constant {N( « const) in this case. These ob- 
servations allow to design easily structures with defined 
spectra and photon-pair emission rates. 



B. Fourth-order interference of a photon pair 

Hong-Ou-Mandel interferometer giving fourth-order 
interference pattern allows to determine entanglement 
time from the profile of fourth-order interference pattern 
obtained by varying a mutual time delay between the 
signal and idler photons. Entanglement time gives a typ- 
ical time window in which an idler photon is detected in 
this specific interferometer provided that its twin signal 
photon has already been detected. Entanglement time 
does not bear any information about the instant of a 
signal-photon detection. Entanglement time defined this 
way does depend on dispersion present in the paths of 
the signal and idler photons in general. However, if the 
same amount of dispersion is present in the signal-photon 
and idler-photon paths entanglement time is not affected 
[131 ■ In case of cw pumping, entanglement time is even 
not sensitive to even orders of material dispersion in the 
paths of two photons from a pair [3l|. On the other hand, 
also sum-frequency generation can be used to define an- 
other typical time constant called correlation time (de- 
fined in optical-field coherence theory) ^3^. Correlation 
time characterizes a true temporal shape of a biphoton 
(as given by the amplitude of a two-photon temporal am- 
plitude) and is sensitive to dispersion. Careful dispersion 
compensation is experimentally required to observe the 
sharpest temporal features as possible. 

If we consider frequency degenerate case, both photons 
in a pair are indistinguishable (are described by identi- 
cal wave-functions) and we observe visibility 100% in a 
Hong-Ou-Mandel interferometer. Different central fre- 
quencies of the signal and idler fields lead to partial dis- 
tinguishability of two photons and, as a consequence, the 
loss of visibility occurs (see Fig.[3K). Moreover, the signal 
and idler photons propagate at different group velocities 
and occur at the output of the crystal with nonzero rela- 
tive mutual delay, that is determined from the shift of the 
dip position in the interference pattern. Also oscillations 
at the difference of the central signal- and idler-field fre- 
quencies occur. Low visibility is usually observed in this 
case because of complex interference of two-photon am- 
plitudes coming from different layers. We note that the 
inclusion of third- and higher-odd-order dispersion terms 
has lead to small oscillations visible in Fig. [3K- Width 
of the coincidence-count interference pattern i?„(r) is 
mainly given by first-order dispersion properties of the 
crystal (for details see, e.g., [231] ) 

Width At of the coincidence-count rate i?„ increases 
with the increasing number Nl of layers concerning uni- 
form poling (see Fig. [3]d), in agreement with decreas- 
ing signal- and idler-field spectral widths. On the other 
hand, width At increases for smaller numbers Nl of lay- 
ers whereas it decreases for larger numbers Nl of layers 
for chirped poling. This behavior can be related to the 
behavior of the signal- and idler-field spectral widths (see 
Fig. [2b): the wider the spectrum the smaller the width 
At. If number Nl of layers is sufficiently large, the width 
At decreases with increasing values of chirping param- 
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eter C that makes spectral widths wider (see Fig. [S]:). 
Several applications (e.g, quantum optical coherence to- 
mography [l^ , determination of ultrashort temporal in- 
tervals [2011) require as short entanglement times as pos- 
sible. Highly chirped crystals are then needed as sources 
of photon pairs giving the widest possible signal- and 
idler-field spectra. 



C. Entropy of entanglement 

Entropy of entanglement E can be used to quantify 
the strength of mutual spectral entanglement between 
the signal and idler photons. Technically, integral equa- 
tions in Eqs. and ifTS]) are discretized and converted 
this way into the problem of matrix diagonalization. Di- 
vision has to be sufficiently close which is indicated by a 
large number of eigenvalues lying close to zero. Eigenval- 
ues of Schmidt's decomposition do not depend on addi- 
tional phase variations introduced by dispersion present 
in paths of the signal and idler photons. This can be 
deduced from the definition of kernels Kg and Ki js^] 
that are proportional to the reduced statistical opera- 
tors of the signal and idler photons, respectively. These 
phase variations only affect phases of the complex base 
functions ipn and (j>n occurring in the decomposition of 
two-photon spectral amplitude A in Eq. (fT3|) . 

Entropy of entanglement E behaves similarly as signal- 
and idler-field spectral widths. The wider the spectra, 
the higher the entropy E. This means, that the higher 
the number Nl of layers in a uniformly poled crystal the 
lower the entropy E. Assuming chirped poling, entropy E 
decreases until a certain number Nlo of layers is reached 
and then increases (see Fig. |4]). Cooperativity parame- 
ter K behaves similarly (see Fig. [4]) as entropy E. We 
note that the wider the signal- and idler-field spectra the 
larger the number of eigenvalues A„ needed in Schmidt's 
decomposition and as a consequence greater values of co- 
operativity parameter K and entropy E occur. 



V. SPATIAL PROPERTIES OF PHOTON PAIRS, 
CORRELATION AREA 



A. Intensity cross-section of the down-converted 
fields 



A typical cross section of the signal-field photon num- 
ber Ns for uniform poling is displayed in Fig. [SK where 
p-polarized signal and idler fields degenerate in frequency 
have been assumed. This means that the signal and idler 
photons propagate as extraordinary waves. The signal 
and idler fields propagate in general as a mixture of TE 
and TM waves at different emission angles, but the dif- 
ference between TE and TM waves is small because of 
small values of radial emission angles and so the pattern 
of cross section is nearly spherically symmetric and oc- 
curs for radial emission angles -ffs lower than 5 deg. The 



b) 





100 200 300 400 500 600 700 800 900 1000 
N, 




c) 

FIG. 3: a) Normalized coincidence-count rate Rn in Hong- 
Ou-Mandel interferometer as a function of relative time de- 
lay T for spectrally degenerate (D) and non-degenerate (N) 
down-conversion in a uniformly poled crystal, b) Width At 
(FWHM) of the interference dip in normalized coincidence- 
count rate Rn as a function of number Nl of layers for uni- 
form (U) and chirped (C) poling, c) Width At as it depends 
on chirping parameter Frequency degeneration is assumed; 
Nl = 1000, C = 1 X 10-" 



highest numbers N of generated photon pairs are found 
in the vicinity oi i)s =0 deg because the structure is 
optimized for collinear interaction. 

Comparison of Figs. [SK and b reveals that chirping of 
poling practically does not change the emission angles. 
It mainly reduces the number N of generated photon 
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FIG. 4: Entropy E and cooperativity parameter K as they 
depend on number Nl of layers. Frequency degeneration as 
well as chirped poling are assumed; = 1 x 10~® fim^^. 

pairs in the middle of the emission pattern, i.e. at the 
emission angles at which the conditions for constructive 
interference of fields from different layers are worse in 
case of chirped poling. As a consequence, a relatively 
wide plateau replaces a typical peak found for uniform 
poling (see Fig. [5b for cuts of the cross-sections). We 
note that the idler twin of a signal photon is emitted into 
a direction given by phase-matching conditions both in 
the transverse plane and along the z axis. 

Signal- and idler-field emission spectra are broad for 
all emission angles as Fig. [6] showing signal-field energy 
spectrum Ss documents. There occurs splitting of the 
signal-field energy spectrum Sg for radial emission angles 
i?s greater than 2 deg. The use of frequency filters thus 
allows to modify the emission-cone cross-section. For ex- 
ample, removal of the central part of the energy spectrum 
Ss leaves the emission-cone cross section in the form of a 
concentric ring. 



B. Correlation area of photon pairs 

Correlation area determines possible emission direc- 
tions of an idler photon [described by radial (di) and 
azimuthal (i/'i) emission angles] which twin, i.e. signal 
photon, has been observed in a given direction defined 
by emission angles -ffs and V's- Because the crystal is not 
infinitely long and also the pump beam is not infinitely 
broad entanglement area has finite width that can be 
characterized by typical widths of the distribution giving 
probability of emission of an idler photon in a given di- 
rection. We denote them as A-i?i and Aipi in the radial 
and azimuthal directions, respectively. Qualitatively, the 
longer the crystal the smaller the correlation area. Also 
the broader the transverse profile of the pump beam the 
smaller the correlation area. 

There occurs splitting of correlation area into two parts 
in radial direction for emission angles that are farther 
from the optimum emission angle -dso for which the struc- 
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FIG. 5: (Color online) Contour plot of the cross-section of 
the signal-field photon number A'^s projected onto hemisphere 
with radial emission angle i9s and azimuthal emission angle ips 
(determines rotation angle around the center of the graph) for 
a) uniform and b) chirped poling, c) Cut of signal-field photon 
number Na along the y axis {ips — deg) traced by radial 
emission angle i^s for uniform (U) and chirped (C) poling. 
Frequency degeneration is assumed; C = lxlO~^ fj,m~^ , Nl = 
1000. 



ture was designed. Considering uniform poling and struc- 
ture optimized for collinear interaction, this splitting is 
clearly visible in Fig. [7^ for radial emission angles iJs 
greater than 0.5 deg. Splitting in radial emission angle 
is accompanied by splitting in frequencies, i.e. frequency 
spectra of each part are separated (compare Figs. [7b and 
[G]) . Chirping of poling leads to blurring of correlation ar- 
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FIG. 6: (Color online) Contour plot of signal-field energy 
spectrum Ss as a function of signal-field emission angle i?3 as- 
suming chirped poling and degenerate spectrum; = 1 x 10~® 
fim-^, Nl = 1000. 



documented in Fig. [7b. This is natural, because 
chirping provides broader spectra of the signal and idler 
fields that allow to fulfill phase-matching conditions in a 
wider area of emission angles. As is evident from Fig.[7h 
chirping partially conceals splitting of correlation area. 
Because width of correlation area is a monotonous func- 
tion of chirping parameter C (see Fig. [HI), photon pairs 
with the required width of correlation area can be gener- 
ated from poled crystals with appropriate chirping. We 
note that width of correlation area is an important pa- 
rameter whenever wave-functions of two photons overlap. 



C. Temperature tuning of the intensity cross- 
section 
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FIG. 7: (Color online) Contour plots showing the cut of cor- 
relation area along the idler-field radial emission angle i9i 
{■di — 'd°^^ + S'&i, is given by optimum phase-matching 
conditions corresponding to a given i?s) as it depends on 
signal-field radial emission angle Frequency-degenerate 
emission of p polarized signal and idler photons along the y 
axis is considered both for a) uniform and b) chirped poling; 
Nl =1000, C = 1 X 10-*^ fj,m-^. 



Properties of the emitted photon pairs in the transverse 
plane can be controlled by temperature of the crystal as 
has been experimentally demonstrated in for chirped 
periodically poled stoichiometric LiTaOs. The change of 
temperature causes changes in values of indices of refrac- 
tion and also widths of layers undergo a certain temper- 
ature shift. As a consequence, temperature can be used 
to tune a sample into resonance where phase-matching 
conditions are perfectly fulfilled. If temperature changes 
from its optimum value, shift in emission directions is ob- 
served (see Fig.[9l). If the crystal is optimized for collinear 
geometry, a spot occurring under optimum conditions is 
replaced by a concentric ring. The greater the declina- 
tion of temperature from its optimum value the greater 
the diameter of this ring. If the crystal is highly chirped, 
this ring is smoothed and a disc is observed. Width of 
the spot as well as width of the ring depend strongly on 
bandwidth of frequency filters placed in front of detec- 
tors. The wider the bandwidth, the wider the spot and 
ring. 
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FIG. 8: Width A&i of correlation area in radial direction 
as a function of for collinear geometry and frequency non- 
degenerate case. The signal and idler photons are p polarized; 
Nl = 1000. 
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investigations of their rich physical properties. 
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FIG. 9: (Color online) Contour plots giving the signal-field 
photon number Ns along the y axis parameterized by radial 
emission angle i9s as it depends on temperature T. The signal 
and idler photons are p polarized. Chirped poling with a) 
= 1 X 10~^ fj.m~^ and b) C = 1 x 10~^ fim~^ is considered; 
Nl = 1000. Spectrally degenerate photons are filtered with 
narrow-band filter 14 THz wide. 



VI. CONCLUSIONS 



A comprehensive study of physical properties of en- 
tangled photon pairs generated in chirped periodically- 
poled nonlinear crystals has been done using a rigor- 
ous model of spontaneous parametric down-conversion. 
Spectral, temporal, and spatial characteristics have been 
addressed. Ultra-wide spectra, high photon fluxes as 
well as entanglement areas can be conveniently controlled 
by chirping parameter. Ultra-wide spectra are respon- 
sible for high values of entanglement of photons in a 
pair that can be useful for quantum-information process- 
ing. The experimentally observed intensity profiles of 
the down-converted fields in the form of discs or rings 
with temperature-dependent parameters have been the- 
oretically explained in this model. It has been revealed 
that complex interference of photon pairs from different 
layers leads to splitting of correlation area that is ac- 
companied by spectral splitting. The presented results 
have shown that chirped periodically-poled crystals are 
not only promising as useful sources of photon pairs for 
many applications but they also remain interesting for 
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